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ABSTRACT 


This paper considers the feasibility of an experiment which 
is intended to measure the density field about a model Helicopter Itotor 
using holographic interferometry. The ei^riment is studied through 
a numerical simulation which we describe after a brief review of the 
measurement technique itself. We then present the data generated by 
this simulation, and assess the prospects both for determining the 
density field from these data, and for actually obtaining such data 
in practice. A few significant problems which may be expected to 
arise are indicated and discussed. 
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INTRODUCTION 


The U.S. Army Aeromechanics Laboratory at NASA-Anes is 
currently planning an experiment to measure the density field 
about a helicopter rotor using holographic interferometry. This 
measurement is to be made in the Laboratory's Anechoic Hover 
Chamber using a 1/7-scale UH-1 rotor. The scale of this test, 
and the character of the flow field to be measured, make this 
probably the most ambitious application of holographic 
interferometry ever contemplated: consequently, it is by no means 
obvious that the application will be successful. Since there is 
a risk involved, some indication that the technique uill work 
uould obviously be desirable before committing substantial 
resources to the project. We at Stanford have therefore 
undertaken an analysis of the feasibility of this experiment. 

Here ue report the results that ue have obtained to date. 

Me begin by reviewing the experimental procedure. 

The first step in making the measurement involves recording a set 
of projections of the flow field. In our case, a projection is a 
profile of integrals of refractive index (which is linearly 
related to density) along rays through the flow field. Thus, the 
value of the projection at point (^.z) and orientation O. as 
shown in figure 1. is 

where n(x.y.z) is the refractive index at (x.y.z). and s is the 
path of the light ray incident at (J.z). The path s is in 
general curved as a result of refraction, although in many 
practical cases its curvature may be neglected. 

Holographic interferometry is used to obtain the set of 
projections as follows. First, a reference hologram of the 
undisturbed field viewed at angle 0 is recorded. Then, a second 
hologram of the disturbed field is made using the same optical 
system. Simultuneous holographic playback of the two produces a 
pattern of interference fringes. The fringe shift at a point 
(^.s) on the pattern is directly proportional to the difference 
between the optical pathlengths to that point through the 
disturbed and undisturbed fields. Thus, if the fringe shift at 
is N. then the change in pathlength is N>s where A is the 
wavelength of the probing light. Since the undisturbed 
refractive index field is known a priori, this interference 
pattern determines the projection of the disturbed field at angle 
0. The same procedure is repeated for a number of viewing angles 
to compose a complete set of data. 
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Uith these data in hand> the object refractive index field 
may be computed using techniques of reconstruction from 
projections. A variety of methods have been formulated for 
performing this calculation, so that reconstruction in oases 
where first, a complete set of projections spanning 180 degrees 
in e is available, and second, the probing rays are straight 
lines, is a matter of routine. If either of these conditions is 
violated, however, the reconstruction problem is considerably 
more difficult. Unfortunately, refraction will not be small in 
our easel we shall discuss the consequences of this presently. 

For the moment, we simply point out that the experiment 
will involve two steps, both of which must be performed 
effectively in order to produce an acceptable result. First, the 
interferometric measurement must make available projection data 
which are sufficiently complete and noise-free to effect 
reconstruction. Second, thr reconstruction method must yield a 
solution having satisfac cry accuracy and resolution. 


THE ASSESSMENT PROCEDURE 

Ue have chosen to assess each of these steps through a 
numerical simulation of the proposed experiment. The input to 
the simulation is a solution for the velocity field about the tip 
of the rotor blade - the region of interest to us - computed by 
Caradonna MI. Ue calculate the corresponding refractive index 
field, and determine its projections by numerical integration. 
These data may then be used as input to a reconstruction routine, 
whose solution may be compared to the original input field. This 
comparison allows us to quantitatively assess the results, which 
is the great advantage of simulation. Thus, we can determine 
precisely the effects of a number of different factors which ue 
may introduce, including noise of various types and magnitudes, 
variations in data sampling rates, different reconstruction 
algorithms, and so forth. 

In particular, ue can expect to obtain the following 


information 

from the simulation: 


1 . 

An impression of the general character of 
project ions . 

the 

2. 

An estimate of the magnitude of the fringe 
be expected. 

shifts to 

3. 

A measure of refraction effects. 


M. 

Minimum sampling rates in ^.z.and « which 
permit acceptable reconstruction. 

will 

5. 

Identification of the most effective reconstruction 
method for the field of interest. 


Ue shall now discuss the simulation and the results obtained thus 
far. 


THE INPUT VELOCITY FIELD 


Let us begin by revieuing the input data. Ue have used 
Caradonna's solution ior a case typical of those to be studied* 
that of zero rotor thrust and a blade tip mach number of 0.9. 

This solution is given as values of the chorduise velocity 
component on a grid of 50 points chorduise by 62 radially by 32 
vertically* spanning the volume 

-5.8 c < K < 10.2 c chorduise* uhere c is the blade chord 

0.75 R < y < 1.75 R radially* uhere R is the blade radius 

0.015 c < z < 22.6 c vertically. 

Most of the significantly disturbed region of the flou field is 
uithin this volume. 

The sensitivity of the interferometric measurement increases 
uith the scale of the disturbed field* so ue must be specific 
about its dimensions in doing the simulation. Ue have used the 
chord and radius of the blade to be tested: 

c = 3 in. =7.62 cm. 

R = *n in. = 1.04 m. 


The corresponding ref motive index field 


For the simulation ue need the refractive index distribution* 
uhich is obtained from the velocity data by first calculating the 
density at each point* and then expressing this as a refractive 
index. The first step involves a form of Bernoulli's equation 
uhich applies for compressible flou that is steady uith respect 
to a rotating (i.e. blade fixed) reference frame. This is 
presented by Caradonna (II: 

<2*- = y Fjzt'. y - 7 


uhere a is the speed of sound* and V is the fluid vector velocity 
measured uith respect to fixed space . The flou is sufficiently 
isentropic that the speed of sound may be eliminated in favour of 
density. Upon doing so ue obtain * 
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The Gladstone-Dale constant K is a iunction oi the refracting 
medium and* ueaki)^* of the wavelength of the probing light. Ue 
intend to use a ruby laser* whose wavelength is 694.3 nm.» in the 
eKperiment; the Gladstone-Dale constant for air at this 
wavelength is 2.2 Sk 10-4 m^/Kg . The reference density and speed 
of sound must also be specified before we can compute the 
refractive index; we have used standard sea level values* namely 


Direct application of these formulae allows us to convert 
from Ue,. the chorduise velocity, to (n-n.r>)* the change in index 
of refraction between the disturbed an undisturbed fields* at 
each point oi Caradonna's solution. Contour plots of several 
horizontal planes through this (n-n^) field are presented in 
figures 2-6* while larger scale views of the region near the 
blade tip are shown in figures 7-11. It is appropriate for us to 
leave the refractive index in this form* rather than to compute 
the disturbed field explicitly, since* as we shall see* each 
inter ferogram may be regarded as a p'ojeotion of the (n-n^) 
field. 


Me should say a few words at this point about the quality of 
these input data. Obviously* if they were known to be correct* 
then we would not be contemplating a difficult experimental 
measurement of them! In fact, while hot-wire data agree well 
with the numerical solution near the blade, they differ 
significantly in the region off the blade tip. The latter region 
is of particular interest in the planned investigation since it 
contains a propagating shock wave. 

However, despite these discrepancies, the fidelity of the 
solution is sufficient for our purposes. For example* while the 
propagating shock is not faithfully reproduced by the numerical 
solution* some of its characteristics are evident in the band 
which may be seen extending off the blade tip in the figures. Ue 
require no more than such general similarity to assess the 
feasibility of the experiment. 


CALCULATION OF FRINGE PATTERNS 


Me shall discuss the refractive index field in more detail 
presently in considering the relationship between the field and 
its projections. For the moment* we shall concentrate upon 
calculation of the projections themselves. 

Me have assumed that refraction is negligible in making 
these calculations. (This is a very desirable situation but we 
cannot justify the assumption solely on the basis of optimism* 
rather, we shall assess its validity after the fact.) Thus we 
have treated the probing rays as straight. Furthermore* we have 
considered them to be parallel to the rotor disc, consistent with 
our plan to use a horizontal beam in the actual experiment. 

In this refrnctionless limit* two substantial simplifications 



emerge. Firsts since rays do not cross horizontal planes through 
the iieldt each plane may be treated independently both in 
calculating the projections and in performing the reconstruction. 
Second, since each ray of probing light traverses the same path 
through both the undisturbed and the object fields, a holographic 
interf erogram will record the projection of the (n-n«) field. 

That is, the fringe shift recorded at a point (%,z) on an 
interf erogram exposed at viewing angle & will be 
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Ue can therefore prrduce a computer-generated interf erogram for 
the viewing angle O by numerically line integrating through 
horizontal planes of our (n-n^) field to a grid of points ('i-.z). 

A set of interf erograms can be composed by repeating these 
calculations for a number of viewing angles. 

Actual interf erograms would of course differ from these computer 
generated versions, since the refractive index field used ai. 
input to our calculations is neither completely accurate nor does 
it encompass the entire region of disturbance - it does not 
include the inner 75% of the blade. But ue have previously 
concluded that the first limitation - inaccuracy of the numerical 
solution - is tolerable, since only approximate information on 
the object field is needed for our feasibility study. For the 
same reason, ue can also do without data on the relatively small 
disturbances over the inner portion of the blade. This is best 
appreciated by examining the fringe patterns themselves, to which 
ue nou turn. 

ORIGINAL PAGE W 

THE CALCULATED FRINGE PATTERNS OP POOR QUALFH' 

Ue have calculated projections for a variety of viewing angles 
between -75 and -i-90 degrees. Figure 12 povidcs a guide to these, 
showing the orientation of each projection with respect to the 
blade. the horizontal coordinate in the projection plane, is 

measured'^ f rom the projected rotor hub: the convention is 
indicated in the figure. Figures 13-22 shou the projections 
themselves. In that ue have plotted these by drawing the loci of 
half-integer fringes (H = tO . 5 , i 1 . 5 , ±2 . 5 , . . . ) and leaving the 
intevening spaces blank, these figures represent so-called 
'infinite-fringe* in ter f erog rams of the flow field. 


All of the projections plotted in this way are quite striking 



and aoma have a oartain aa&thetio appeal. But perhaps the most 
impressive of all it the chorduise projection (figure 13). Fori 
when plotted as an infinite fringe interference pattern, it is 
completelp blank. Nowhere, apparently, do the optical 
pathlengths through the undisturbed and object fields differ by 
as much as one-half wavelength - 350 nanometers. An actual 
interf erogram would therefore consist only of diffuse gray bands 
and would have very little contrast. 

This result may be understood by considering the plots 
of the (n-n^) fields in figures 2-11. Note that the region of 
significant disturbance is thin in the chorduise direction, 
extending little beyond the leading and trailing edges of the 
blade. Thus, disturbances over only a short segment of any line 
through the field contribute to the projection integral, 
noreover. the sign of the disturbance oscillates over this short 
segment from positive ahead of the blade, to negative in the 
region of reduced pres'ure above it. and back to positive behind 
the trailing edge. Therefore, disturbances over the middle of 
the line segment tend to cancel those at either end. and as a 
result the net line integral is very small. 

At the other extreme, we have the radial (90*) projection. 

The disturbed region is long in this direction, and most of the 
disturbance along any given line through the field is of the same 
sign. Thus large fringe shifts are produced - a bright fringe of 
order -*-25 lies just ahead of and above the projected leading edge 
(5=0). while one of order -25 lies slightly above the 1/4-chord 
of the blade. Integration over the inner IS*/, of the blade, which 
has been excluded from our calculations, would tend to increase 
these magnitudes. 

There is. then, a large gradient in fringe shift across the 
radial projection. This is troublesome. An immedeately apparent 
problem is that individual fringes become indistinct in such a 
strong gradient; thus in order to accurately determine the fringe 
order over this rather important region of the projection one 
would somehow hp.vc to either reduce the sensitivity of the 
measurement or improve the clarity of the fringe patterns. A more 
fundamental problem is that such closely spaced fringes indicate 
that refraction is most certainly not negligible, as we can 
demonstrate by the following calculation. Vest 12] suggests the 
criterion that if the quantity 



exceeds about 0.3. then refraction should be taken into account. 
Here d^. is the fringe spacing, and L is the distance between 
imaging lenses in the optical system. In our case, we have 

C. . ^ /‘C // . 


Therefore, according to this criterion, refraction will be 
significant if 
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Fortunately, uhile this condition is violated by about 1 
order oi nc'initude in the radial proieotion. it is at least 
Marginally satisfied at lower angles. Consider, for ONaMpie. the 
projection at +75 degrees. Rays traversing the field at this 
angle are nearly tangent to the band of disturbance extending off 
the blade. He would expect this fairly long transit through a 
band of significant variation in the refractive index to produce 
a large fringe shift, and indeed this is the case. The fringe 
order varies from *-17 to +1.5 across the projection. (These 
values would not change very much if the field over the inner 75X 
of the disc were included, since rays within the band off the 
blade have only short segments within the disturbed part of that 
field.) The minimum fringe spacing in this projection is 3~b mm., 
which may be considered marginally satisfactory with 
respect to the refraction criterion. 

As the projection angle is reduced further from chordwise. 
the fringe spacing continues to increase, and the magnitude of 
the fringe shift correspondingly decreases. This is a 
consequence of both shortened pathiengths through the disturbed 
region and oscillations in the sign of (n-nj.) along these paths. 
Similar behavior is evident for angles between -75 degrees and 0| 
however, the fringe pattern for any given negative projection 
angle is less rich than that for its positive counterpart, since 
rays inclined at negative angles travel across rather than along 
the disturbed band off the blade tip. An important feature here 
is the pronounced variation in the nature of the fringe patterns 
with projection angle. This characteristic of the (n-no) field 
bears significantly upon the field of view required for 
reconstruction, as we shall see. 


THE FEASIBILITY OF RECONSTRUCTION 


Having calculated and examined the projections of the 
refractive index field, we can now address the following 
question: can this field be reconstructed with satisfactory 

accuracy and resolution using input data from these projections? 
Although we have not yet attempted a reconstruction, we can 
indicate some of the factors which will affect its feasibility. 

The theory of reconstruction from projections affords an 
analytical solution provided that two conditions are satisfied: 

1. Complete knowledge of the projection. over 

a 180 degree range of viewing angles is available. 

2 . Refraction is negligible. 

In practice, the analytical solution is adapted to accept disrete 
data. Reconstruction algoriths differ in their input 
requirements - some can accept input points spaced non-uni f orm 1 y 
in N and 0. for example - but a guide is provided by the sampling 
theorem. This theorem indicates that input data points should be 
separated by no more than twice the maximum spatial frequency of 
the disturbance field in the direction perpendicular to the 
projection! otherwise aliasing will produce artifacts in the 
reconstruction. This is a lower limit, however, since redundant 


data ara 


aluaya raquirad to 


auppraaa noiaa 


Nout in aampiinp Iringa pattarna* ona uannot dataraina 
praoisaly tha iringa ordar at aach oi a aat oi uniic uly apaoad 
pointa. Ona can only dataraina uhathar a aaleotad point liaa within 
aithar a bright ( integar-*ordar ) or a dark (hal f-intagar ordar) 
iringa* and tha rasolution oi a diraot aaaauraaant ia tharaiora 
liaitad to 2 1/H iringa ahiit. For anaapla* ona would aaaign tha 
value zero to aach point on a iringa pattern auch aa that ahown 
ior tha chordwiaa viewing direction* although tha actual value 
variaa between ^0.2. 


There ara techniquea whereby aharper raaolution can be 
obtained. Tha ainpleat way would be by graphical interpolation 
between the oentrelinea oi bright and dark iringeai thia 
technique would be applicable to the richer iringe patterna 
obtained at higher projection anglea. Another method would 
involve the uae oi so-called ' i ini te-i r inge* interierograna . 
Thua*one intr?»ducea an artiiicial pathlength gradient acroaa the 
projection between holographic exposures oi the reierenca and 
object iields. Then ii the object iield contained no 
disturbance* holographic playback oi the two exposures would 
produce an interierence pattern consisting oi equally spaced 
parallel iringes. The actual disturbance iield will displace 
these iringes* and the pathlength data is extracted by measuring 
their displacement. Again* the sampling rate ia determined by 
the iringe spacing* but by using this technique that spacing can 
be reduced arbitrarily - at least to the point where adjacent 
iringes become indistinct. 

While it is possible* then* to ret Ive iringe shiits to 
higher accuracy than t. 1/2* this level oi accuracy should be 
quite acceptable ior input to a reconstruction routine. Moreover* 
higher accuracy would perhaps be unjustiiied in view oi the 
noise level to be expected. 


Field oi view 

Oi course* the iield oi interest to us contains shock 
waves* which are very high irequency structures. Thereiore* we 
can anticipate that quite high sampling rates in both and 0 
will be required in projections near tangent to the shocks - 
those within about 20 degrees oi spanwise - to obtain an 
alias-iree reconstruction having reasonable spatial resolution. 

This suggests that the present apparatus will place a very 
severe limitation on th** quality oi the reconstruction* ior it 
has a restricted iield oi view which will prevent recording oi 
projections near spanwise. Th i v restriction can be 
appreciated by considering iiguie 23. The probing light passes 
across the test chamber* irom one side oi the rotor disc to the 
other. li only one blade's disturbance iield is to be included 
in the reconstruction process (the alternative would obviously 
increase the experimental and computational burden enormously)* 
then each ray which traverses one blade's iield must not 
intersect that oi the other blade. Since the region oi 
signiiicant disturbance extends outward irom about haliway along 
the blade* and about 1.S chords out irom the leading and trailing 
edges (the iormer oi which is itseli about 1/4 chord ahead oi the 





dianeter oi the rotor disc) this restriction excludes angles 
within about 20 degrees oi spanuise. 

While this limitation oi the iield oi view to about tOO degrees 
can be tolerated ii the object field is 'roundish* (ii not 
completely axisymmetric) - provided that clean projection data is 
available and that care is used in the reconstruction - ior our 
iield. whose character varies drastically with viewing angle, a 
iull 180 degree iield oi view will likely be required. This 
could only be obtained with an experimental arrangement diiierent 
irom that presently contemplated. Three possible schemes are as 
iol lows : 

First, we might place a mirror on the rotor hub to 
reflect the beam away irom the second blade. However, the 
problems of mounting a mirror of the necessary size on the hub 
while meeting stringent requirements ior optical quality and 
stability are formidable; therefore this scheme would probably not 
be practicable. 

Second, we might rearrange the optical system so that the 
probing beam was parallel to a vertical plane rather than the 
rotor disc. Ue could measure only that portion of the 
disturbance field outside the disc using this arrangement, since 
blockage by an opaque blade, or refraction by a transparent 
blade, would preclude recording of projections within the 
diameter of the disc. This scheme would have additional 
disadvantages in that it would mandate extensive modifications to 
the test chamber, and the entire optical system would have to be 
moved to record each new projection. The results would probably 
not justify the effort involved. 

A third possibility would be to use a s i ng 1 e-b 1 aded rotor. 

Such a rotor would place quite a strain on the dynamic components 
of the rotor system, and might require some modifications to the 
affected hardware. While we would prefer to avoid the consequent 
effort and expense of implementing this scheme, no practical 
alternative which would satisfy the field of view requirement is 
evident. 


Ref raction 


Another significant problem which we have mentioned earlier 
is refraction. Since refraction is not negligible in 
projections within about 25 degrees of spanwise. the 
r ef ract i on 1 ess limit cannot be assumed in performing the 
reconstruction. T!>e question is rather whether a first-order 
correction will take this effect into account satisfactorily, or. 
instead, a full-blown iterative reconstruction will be necessary. 
A discouraging sign is that both the horizontal and the vertical 
fringe spacings in the radial projection violate the criterion 
for negligible refraction. This implies that horizontal planes 
cannot be treated independently in the reconstruction, and that a 
three-dimensional scheme will be necessary. Reconstruction in 
such circumstances is by no means a closed problem and a 
satisfactory solution consequently cannot be assured; in any 
case, the difficulties of developing a reconstruction scheme and 
of implementing it would be substantial. 


NOISE IN THE EXPERIMENTAL DATA 


For the monent let us assume that a practical reconstruction 
scheme can be developed> and consider the question of whether 
sufficiently noise-free input data for such a routine can 
actually be obtained in practice. Actually, noise may 
well make these data rather difficult to measure. Me can 
illustrate this by a simple calculation. 


Variations in temperature between recordings of the 
reference and object exposures will be one source of random 
noise. Let us estimate the temperature stability which 
would have to be maintained in order to measure the same 
fringe shifts that ue have calculated in our simulation. 
Assume that, between exposures of the reference and object 
fields, the temperature changes from To to T. Then the 
fringe shift produced by this temperature change alone is 
given by 
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where L is the physical pathlength in the optical system. dn/dT 
may be found using the equation of state and the Gladstone-Dale 
formula. Me have that 


^ -7 

o7', 


Of 


V" 


From the equation of state. 


Of _ 

c<r 


77 r 


and . 

from the 

Gladstone-Dale 

formula . 


/ 

- A A 

/C 

A. 

Thus . 





c/n 
r<:' 1 " 

/ 

# 


and 


' ■■ X - > 

/ , ^ 


71 a A- 



A . 


for L in metres and dT in degrees Kelvin. In our case. L is 
about 20 m. If uc require that temperature variations produce 
fringe shifts of no more than 1/2. tiien we must maintain constant 
temperature to within 


t" A d>'yyC2.e:)) 

NoUf a unilorn tenp'drature ohanga baiuaen aaposuras of tha 
raierenco and object fielus of this magnitude* or any other for 
that matter* could be detected easily and then removed from 
the measurement. However* random variations in temperature 
across the field cannot be dealt with so easily. While random 
variations would have to be of larger amplitude to produce a 
1/2-lringe shift noice level in the projections* such 
variations are not inconceivable since air is continuously 
pumped from outdoors into the Anechoic Room* and around inside 
the facility* while the ewperiment is in progress. This effect 
must therefore be of some concern. 

Nor is this the only source of noise. Besides additional 
sources of random noise* there is the problem of unsteadiness in 
the flow field. All projections of such a large field cannot be 
measured simultaneously; in fact* since the region of significant 
disturbance is wider in some directions than the aperture of the 
optical system* even single projections must be composites. In 
order to compose a projection the optics must be moved with 
respect to the rotor* while to assemble a set of projections 
exposures must be made with the rotor at various different 
azimuths. Consequently* since there is both a periodic variation 
in the disturbance with azimuth due to the non-axisymmet ry of the 
Anechoic Room* and random unsteadiness* different parts of some 
projections* and projections taken at different angles* will in 
fact be projections of somewhat different flow fields. In this 
circumstance the best that one could hope to obtain fiom the 
reconstruction would be a rather blurred image of the flow field. 
High frequency components* particularly shock waves* would tend 
to be washed out of the image* since these features are of 
primary interest* the utility of the results would therefore be 
much reduced. 


SUMMARY 


Through our simulation and analysis of the planned experiment 
we have identified a number of factors which may be expected to 
complicate the holographic measurement or limit the quality 
of the density field reconstruction. These include: 


1 . 

The restricted 

field of view of 

the present apparatus. 

2. 

Refraction in 

projections near 

spanwise. 

3. 

Random noise. 



4. 

Unsteadiness in the flow field. 



Me intend to do a number of reconstructions from our simulated 
projections which will provide further information on the effects 
of these factors. The cases to be considered are as follows: 


1. Reconstruction with a full 180 degrees of noise-free 
projection data. 



2. Reconstruction with the field of view limited to 
mo degrees. 

3. Reconstruction with e 180 degree field of view and 
artificially generated random noise superimposed 
upon the projection data. 

R. Reconstruction with both superimposed noise and a 
restricted field of view. 

Comparison of the reconstructions calculated in the last three 
cases with that of the first case uill allow quantitative 
definition of the limitations imposed by noise and a restricted 
field of view. Ue shall present these results in a later report 
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Figure 14. PROJECTION OF THE (n - n ) FIELD AT +30® FROM CH0RDWI5E 
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FlKure 18. RADIAL PROJECTION OF THE (n - n ) FIELD 
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